1 



General temporal instability criteria for 
stably stratified inviscid fiow 

By LIANG SUN^^f 

^School of Earth and Space Sciences, University of Science and Technology of China, Hefei, 

230026, China. 

■^Dept. of Modern Mechanics, University of Science and Technology of China, Hefei, 230026, 

China. 

(Received 20 April 2010, and in revised form 19 October 2011) 

The temporal instability of stably stratified flow was investigated by analyzing the Taylor- 
Goldstein equation theoretically. According to this analysis, the stable stratification N'^ ^ 
has a destabilization mechanism, and the flow instability is due to the competition of 
the kinetic energy with the potential energy, which is dominated by the total Froude 
number Fr^. Globally, Fr^ ^ 1 implies that the total kinetic energy is smaller than 
the total potential energy. So the potential energy might transfer to the kinetic energy 
after being disturbed, and the flow becomes unstable. On the other hand, when the 
potential energy is smaller than the kinetic energy {Fr^ > 1), the flow is stable because 
no potential energy could transfer to the kinetic energy. The flow is more stable with the 
velocity profile U' /U'" > than that with U' /U'" < 0. Besides, the unstable perturbation 
must be long- wave scale. Locally, the flow is unstable as the gradient Richardson number 
Ri > 1/4. These results extend the Rayleigh's, Fj0rtoft's, Sun's and Arnol'd's criteria 
for the inviscid homogenous fluid, but they contradict the well-known Miles-Howard 
theorem. It is argued here that the transform F = (j>/{U — c)'^ is not suitable for temporal 
stability problem, and that it will lead to contradictions with the results derived from 
the Taylor- Goldstein equation. However, such transform might be useful for the study of 
the Orr-Sommerfeld equation in viscous flows. 



1. Introduction 

The instability of the stably stratified shear flow is one of main problems in fluid 
dynamics, astrophysical fluid dynamics, oceanography, meteorology, etc. Although both 
pure shear instability without stratiflcation and statical stratification instability without 
shear have been well studied, the instability of the stably stratified shear flow is still 
mystery. 

On the one hand, the shear instability i s known as the instability of vortici ty max- 
imum, after a l ong way of investigations (jRavleigh II88OI : iFjartoftl I1950I : lArno ld 19651 
SunI It is recognized that the resonant waves wit h special v elocity of the 

concentrated vortex interact with flow for the shear instability ( Sun 20081) . Other veloc- 
ity profil e s are stable in homogeneous fluid without stratiflcation. On the other hand, 
Ravleig proved out that buoyancy is a stabilizing effect in the statical case. 

Thus, it is naturally believed that the stable stratiflcation do favor the stability (see, e.g. 
Taylor 1931 ; Chandrasekhar 196lh , which finally results in the well known Miles-Howard 
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theorem (iMilesl I 196ll Il963t iHowardI 196ll ). Accordmg to this theorem, the flow is sta- 
ble to perturbations when the Richardson number Ri (ratio of stratification to shear) 
exceeds a critical value Ric = 1/4 everywhere. In three-d imensional stratified fiow, the 
corresponding criterion is Ric = 1 ( Abarbanel et al. 19841) . 

However, the stabilizat ion effect of buoyancy is a illusion. In a less known paper, 
Howard fc Maslov3 ( 19731 ) had shown with several special examples that stratification 



effects can be destabilizing due to the vorticity generated by non-homogeneity, and the 
instability depends on the det ails of the veloci t y and density profiles. One insta bility is 
called as Holmboe instability (jHolmbod Il962t lOrtiz et al\ 12002 ; lAlexakislliooi ). Then 
iHoward fc Maslovd (fl97i stated three main points from the examples without any fur- 
ther proof, (a) Stratification may shift the band of unstable wave numbers so that some 
which are stable at homogeneous cases become unstable, (b) Conditions ensuring stabil- 
ity in homogeneous fiow (such as the absence of a vorticity maximum) do not necessarily 
carry over to the stratified case, so that 'static stability' can destabilize, (c) New physical 
mechanisms brought in by the stratification may lead to instability in the form of a pair 
of growing and propagating waves where in the homogeneous case one had a stationary 



wave. 



Recall the points bv iHoward fc Maslovd ( 19731). and tha t there is a big gap between 



Rayleigh's criterion and Miles-Howard' criterion. lYihl ( 1980f ) even wrote "Miles' criterion 
for stability is not the nature generalization of Rayleigh's well-known sufficient condition 
for the stability of a homogeneous fiuid in shear flow" . The mystery of the instability is 
still cover for us. 

Following the frame work of Sun ( 2007 . 20081) . this study is an attempt to clear the 
confusion in theories. We find that the flow instability is due to the competition of the 
kinetic energy with the potential energy, which is dominated by the total Froude number 
Frf . And the unexpected assumption in Miles- Howard theorem leads the contradiction 
to other theories. 



2. General Unstable Theorem For Stratified Flow 

2.1. Taylor- Goldstein Equation 
The Taylor- Goldstein equation for the stratified inviscid flo w is employed (jHowar ll96ll 



Yih 1980t Baines fc Mitsuderal 1994t Criminale et al. 2003), which is the vorticity equa- 
tion of the disturbance ( Drazin fc ReidI 20041) . Considering the flow with velocity profile 
U{y) and the density field p(y), and the corresponding stability parameter N (the Brunt- 
Vaisala frequency), 

= -gp'/p, (2.1) 

where g is the acceleration of gravity, the single prime ' denotes d/dy, and N"^ > denotes 
a stable stratification. The vorticity is conserved along pathlines. The streamfunction 
perturbation cj) satisfies 



U" 



U 



(2.2) 



where k is the real wavenumber and c = -\- ici is the complex phase speed and double 
prime " denotes jdy^ . For k is real, the problem is called temporal stability problem. 
The real part of complex phase speed is the wave phase speed, and = ka is the 
growth rate of the wave. This equation is subject to homogeneous boundary conditions 



(j) = at y = a,b. 



(2.3) 
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It is obvious that the criterion for stabihty is uJi — (ci = 0), for that the complex 
conjugate quantities and c* are also physical solutions of Eq. ()2.2p and Eq. p.3|) . 

Multiplying Eq. ()2.2|) by the complex conjugate </>* and integrating over the domain 
o- ^ y ^ b, we get the following equations 



U"{U - Cr) 

|C/-cP 



\<P\']dy= I "'' N^Wdy. (2.4) 



and 



U" 



2{U ~Cr)N^' 



\U- 



M'dy^O. 



\U-c\^ \U~c\^ ^^'^^ 

In the case of N"^ = 0, iRavleig used Eq. ()2.5p to prove that a necessary condition 

for inviscid instability is U"{ys) = 0, where Vs is t he inflection point and Us = U{ys) 
is the velocity at j/s- Using Ea. (l2.5p . ISvngd (|1933l ) also pointed out that a necessary 
condition for instability is that U" — should change sign. But such condition 

is useless as there are two unknown parameters Cr and c^. 

As a first step in our investigation, we need to estimate the ratio of \(j)'\'^dy to 

\(j)\'^dy. This is known as the Poincare's problem: 



WVdy 



\c^?dy, 



(2.6) 



where the eigenvalue /i is positively definite fo r any 4) # 0. The smallest eigenvalue value, 
namely /ii, can be estimated as /ii > (^^)^ (|Mu et a/.lll994t ISunll200l1) . 



2.2. General Instability Theorem 

In departure from previous investigations, we shall investigate the stability of the flow by 
using Ea. (l2.4p and Eq. (|2.6l) . As /i is estimated with boundary so the criterion is global. 
We will also adapt a different methodology. If the velocity profile is unstable (c; ^ 0), 
then the equations with the hypothesis of = should result in contradictions in some 
cases. Following this, a sufficient condition for instability can be obtained. 
Firstly, substituting Ea. (|2.6p into Ea. (|2.4p . we have 



9{y) 



U- 



'■dy = - 



Kv) 

\U-c\- 



''dy, 



where 



2iV2 



(2.7) 

(2.8a) 
(2.86) 

It is noted that g{y) > for ^ 0. Then cf > if h{y) ^ throughout the domain 
o- ^ y ^ b for a proper Cr and k. Obviously, h(jj) is a monotone function of k: the smaller 
k is, the smaller h{y) is. When fc = 0, h{y) has the smallest value. 



giy) = t^ + k 

h{y) ^{fi + k^){U ~ Crf + U"{U - Cr) - 



\u- 



h{y) = N 



2AU-Crf 



U"{U - Cr) 



If we define shear, parallel and Rossby Froude numbers Frt, Frs and Frr as 



Fr^ = Fri 



Frl,Frl=^^ ^''^ 



p2_ U"{U-Cr) 



(2.9) 



(2.10) 
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where the shear Froude number Fts is a dimensionless ratio of kinetic energy to po tentia l 



energy. As U" plays the same role of /3 effect in the Rossby wave (.Smi 2006 al, 120071 ) , 
the shear Froude number Frj. is a dimensionless ratio of Rossby wave kinetic energy to 
potential energy. Then h{y) ^ equals to Frj ^ 1. Thus a general theorem for instability 
can be obtained from the above notations. 

Theorem 1: If velocity U and stable stratification N"^ satisfy h{y) ^ or Fr^ ^ 1 
throughout the domain for a certain c^, the flow is unstable with a q > 0. 

Physically, Frf ^ 1 implies that the total kinetic energy is smaller than the total 
potential energy. So the potential energy might transfer to the kinetic energy after being 
disturbed, and the fiow becomes unstable. On the other hand, when the potential energy 
is smaller than the kinetic energy {Frf > 1), the flow is stable because no potential 
energy could transfer to the kinetic energy. 

Mathematically, we need derive some useful formula for applications, since there is still 
unknown Cr in above equations. To this aim, we rewrite Ea. (l2.9p as 

Till jjia 

Ky) =fl{U+—- Crf - {N^ + —). (2.11) 

Assume that the minimum and maximum value of within a ^ y ^ bis respectively 

rrii and rria- It is from Eg. (12. Ill) that ^ ^ TOq for the smallest value of h(jj). Thus 
a general theorem for instability can be obtained from the above notations. 

Theorem 2: If velocity U and stable stratification iV^ satisfy h(y) ^0 throughout the 
domain for a certain rrii ^ c,. ^ TOq, there must be a Ci > and the flow is unstable. 

It is from Eq. ipiT)) that h{y) ^ requires /x(C/ + ^ - c^)^ less than {N'^ + ^). 
The bigger is, the smaller h{y) is. So the stable stratification has a destabilization 
mechanism in shear fiow. This conclusion is new as former theoretic studies always took 
the static stable stratification as the stable effects for shear fiows. 

According to Ea. (|2.1ip . the bigger — is, the more stable the fiow is. It is obvious 
that uia - mj is bigger for U'/U'" > than that for U'/U'" < 0. So the flow is more 
stable with the velocity proflle U'/U'" > 0. 

Although Theorem 1 gives a sufficient unstable condition for instability, the compli- 
cated expression makes it difficult for application. In the following section we will derive 
simple and useful criteria. 



3. Criteria For Flow Instability 

3.1. Inviscid Flow 

The simplest fiow is the inviscid shear fiow with N'^ — 0. The sufficient condition for 
instability is h{y) ^ 0. To find such condition, we rewrite h{y) in Eq. (|2.11| ) as 

h{y) = {Ul - Cr)iU2 - Cr) (3.1) 

where Ui = U and U2 — U + U" / fi. Then there may be three cases. Two of them have 
Ul intersecting with U2 at U" = (FigH]). The first case is that U"/{U - U^) > 0; 
thus, h{y) > always holds at = Us as shown in FiglT^. The second case is that 
U" /{U — Us) < 0; thus, h{y) < always holds in the whole domain, as shown in FiglT]D. 
In this case, the flow might be unstable. 

The sufficient condition for instability can be found from Ea. (l3.ip as shown in FiglTjD. 
Given Cr = Us, Eg. p.ip becomes 
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Figure 1. The value of h(y) under the condition U'^ = 0: (a) for i7"/((7 - Vs) > 0, (b) for 

U"/{U - Us) < 0. 
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Figure 2. The value of h{y) for U" / in case 3 and case 4. 



If 



u" 



f^jj_u ^ < — M is alvi^ays satisfied, h{y) < holds within the domain. 

Corollary 1.1: If the velocity profile satisfies J^jj < — /i within the domain, the flow 
is unsta ble. 

Since ISunI ( 2007 ) obtained a sufiicient condition for stability, i.e. iFtt - > ~^ within 
the domain. The above condition for instability is nearly marginal ()Sun 2008 ). 

The last case is that U" ^ throughout the domain; thus, h{y) > always exists 
somewhere within the domain, as shown in Figl2^. 

3.2. Stably Stratified Flow 

If the static stratiflcation is stable (iV^ > 0), then g{y) is positive. The flow is unstable 
if h{y) is negatively deflned within a^y^b at k = 0. We rewrite h{y) as 

h{y) ^ fJ-{Ul ~ Cr){U2 - Cr) 

= ti[U + ^{U" - VU"^ + 4.^iN^ ) - cr] (3 3^ 

x[t/ + ^(L/" + v/C/"2 + 4/iiV2 ) _ c^]^ 

The value of h{y) can be classifled into 4 cases. The flrst and the second ones ([/" = 
and iV^ = a,t y = ys) are similar to discussed above and shown in Fi^T^ and Fi^TJj. For 
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such cases, we have a sufficient condition for instability. 

U"{U - Us) - 



{U - Usf 



< -M- (3.4) 

This can be derived directly from Eq. (|2.7p . similar to Corollary 1.1. The first sufhcient 
condition for instability is due to the shear instability, and the unstable criterion is 
Eq.dSH). 

Corollary 1.2: If the velocity profile satisfies ^ ^ within the domain, 

the flow is unstable. 

The third case ([/" ^ 0) is also similar to the case in Fi^2^, and the flow is stable. 
The last one is unstable flow shown in FiglllD, where U" ^ and h{y) < throughout. 
In the last case, the maximum of t/i must be smaller than the minimum of U2 so that 
a proper Cr within the Ui and U2 could be used for the unstable waves. Although the 
exact criterion can not be obtained as the required maximum and minimum can not be 
explicitly given, the approach is very straightforward. 

Nevertheless, we can also obtain some approximate criterion for the fourth case. It 
is from Eq. (|2.1ip that h{y) ^ if the minimax of /i(J7 + ^ — Cr)^ is less than the 

minimum of (7V^ + ^^). As the minimax value of ^{U + ^ — c^)^ is jfi{ma — mi)"^ 
when Cr — {nia + TOi)/2, we obtained a new criterion according to Ea. (|2.10p . 

Thus a sufficient (but not necessary) condition for h{y) < is that the following equation 
holds for a ^ y ^ 6. 

Fri s; 1. (3.6) 

From the above corollaries, the flow might be unstable if the static stable stratification 
is strong enough. The stably stratification destabilize the flow, which is a new unstabl e 



mechanism. The above corollary contra dicts the previous results (lAbarbanel et aLlll984r). 



but it agrees well with t he recent theor y ( Friedlandeijl2001 ) , experiments ( Zilitinkevich et al. 



2008f) and simulations ( Alexaki^l2009f) . Again, we point out here that the flow is unstable 
due to potential energy transfer to kinetic energy under the condition of Frf ^ 1. 

This conclusion is new because it is quite different from previous theorems in which 
the static stable stratification plays the role as a stabilizing factor for shear fiows. 



4. Discussion 

4.1. Necessary Instability Criterion 

In the above investigation, it was found that stable stratification is a destabilization 
mechanism for the fiow. Such finding is not surprising if one notes the terms in Eq. ()2.2p . 
Mathematically the sum of terms in square brackets should be negative for the wave 
solution. Thus both > and N'^ < do favor this condition. This is why the 

unstable solutions always occur at < in shear flow. And > here might lead 
to cf > 0. Physically, the perturbation waves are truncated in the neutral stratified flow. 
But the stable stratification allows wide range of waves in the p erturbation. Such waves 
might interact with each other like what was illustrated in [SunI ([2OO8I . 

As Theorem 1 is the only sufficient condition, it is hypothesized that the criterion is 
not only the sufficient but also the necessary condition for instability in stably stratified 
flow. This hypothesis might be criticized in that the flow might be unstable (cf > 0) 



General instability criterion 



7 



if h(y) changes sign within the interval (Fi^2^), where a proper chosen (j) would let the 
right hand of Eq. (|2.7p become negative. 

However, this criticism is not valid for the case in Fi^5^. It is from the well-known 
criteria (e.g. Rayleigh's inflexion point theorem) that the proper chosen (/) always let 
the right hand of Ea. (l2.7p vanish. It seems that the flow tends to be stable, or the 
perturbations have a prior policy to let a — 0. The flow become unstable unless any choice 
of (j) would let the right hand of Eq. ()2.7p be negative. In this situation, we hypothesize 
that Theorem 1 fully solves the stability problem. 



4.2. Long-wave Instability 

In inviscid shear flows, it has been recognized that very short-wave perturbations are 
dynamically stable under neutr al stratification, and the dynamic instability is due to the 
larger wavelengths (|Sunll2006a) . It should be noted that Rayleigh's case i s reduced to the 
Kelvein-Helmholtz vortex sh eet model under the long-wave limit fc <C 1 (jHuerre fc Rossi 
1998t ICriminale et a/.l [20031) . We have shown that this can be extended to shear flows, 
and that the growth rate Wi, is proportional to JJii ( Sun 20066 . 20081) . 



Such conclusion can be simply generated to the stratified shear flows, which can be seen 
from Ea. (|2.8p . If k is larger than a critical value kc, the sufficient condition in Theorem 
1 can not be satisfied and the flow is stable. For shortwave (fc ^ 1), h{y) is always larger 
than that for long -wave k <C 1 . The long -wave instabilit y in the stratified shear flow 
was also noted by iMiled (Il96ll Il963l) and E owar 

3 (|l96ll) . who showed a likelihood of 
Ci — > at fc — » oo. The long - wave instability theory can explain the results in numerical 
simulations ( Alexakisll2009f) . where the unstable perturbations are long- wave. 



4.3. Local Criterion 

In the above investigations, an parameter /i is used, which represents the ratio of two 
integrations with boundaries. So the criteria are global. On the other hand, we can also 
investigate the local balance without boundary conditions. For example, consider the 
flow within a thick layer —S ^ y ^ 6. The velocity is U{y) = Uq + U'y, and the kinetic 
energy is (C/ — c^)^. The stratification is N'^, and the potential energy is N'^cP, where 
d — 26 is the thickness of the layer. The Froude number is Fr^ = (U'^S^)/ (N^cP) for 
Or — Uo- The instability criterion in Eq. p.6p becomes 



J7'2 I 

R^-^>1- (4-1) 

If local gradient Richardson number exceeds 1/4, the local disturbances is unstable. 
However, the flow might be stable as the globe total Froude number Fr^ > 1. This 
criterion is opposite to Miles-Howard theorem, we will show why Miles-Howard' theorem 
is not correct from their derivations. 



4.4. Relations to Other Theories 
In th e inviscid shear flow, t he linear theories, e.g., Rayleigh-Kuo cirterion ([Criminale et 



20031 ). Fj0rtoft criterion (iFiortoftI [l950l) and Sun 's criterion (jSunI l2007t ). are equal to 
Arnol'd's nonlinear stability cri terion jArnoldj 19651). Arnol'd's first stabilit y theorem cor- 
responds to Fj0rtoft's criterion ( Drazin fc Reidll2004 : Criminale et a l!^200^ . and Arnol'd's 
second nonlinear theorem corresponds to Sun's criterion ( SunI 2007 . 2008). It is obvious 
that the present theory, especially Corollary 1.1 is a natural generalization of inviscid 
theories. 

In the stratified fiow, iMilesI ([l96ll,[i963) and lHowardI (|l96l[ ) applied a transform F = 
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(j)/ {U — c)" to Eq. ()2.2p . which aUows different kind of perturbations. Thus n = 1/2 gives 
Miles's theory and n — \ gives Howard's semicircle theorem. 

Considering that n — 1 and = (|Howardlll96l . Eq.(l231) becomes 



{\F'\'' + e\F\^)[{U-Crf-4]dy = G. (4.2) 

It is from Eq. (|4.2l) that all the inviscid flows (no mater what the velocity profile U{y) 
is) must be temporal unstable (fc is real). This contradicts the criteria (both linear and 
nonlinear ones) for inviscid shear flow. So the wavenumber k in Eq. (j4.2l) should be com- 
plex k = kr + iki. Besides, from Eq. (|2.11[ ). Eq. p.3p and Figl^jD, the unstable Cr might be 
either within the value of U or beyond the value of U. This also contradicts the Howard's 
semicircle theorem for the stratified flow. It implies that the transform F is not suitable 
for temporal stability problem. 

Taking n — 1/2, Howard extracted a new equation from Taylor- Goldstein equation, 

[{U - c)F'r - [k^{U -c) + ^ + i\u" N')/iU c)]F = (4.3) 

After multiplying above equation by the complex conjugate of F and integrating over 
the flow regime, then the imaginary part of the expression is 

1 

|Ff + [/fc2|F|2 + (-[/'2 - N^)\F\y\U - c|2 = (4.4) 

Miles- Howard theorem concludes that if c,; 7^ 0, then Ri < j for instability. 

However, the transform F — 4>l \/lJ — c requires a complex function F , even though 
both (j) and c are real. In that \/U — c might be complex somewhere as U — Cr < 0. 
Consequently, the wave number k in Eq. (|4.3|) is a complex number but no longer a 
real number as that assumed in Taylor- Goldstein equation. The complex wavenumber 
k leads to spatial stability problem but temporal stability problem investigated in this 
study. The assumption of = with ki ^ Q implies the flow is unstable with 7^ 0. 
However, Howard ignored this in his derivations. That is why Miles-Howard theorem 
leads contradictions to the present studies. 

Although the transform F = 0/([/ — c)" leads some contradictions with Rayleigh crite- 
rion and present results, it might be useful for the viscous flows. In these flows, the spatial 
but temporal stability problem is dominated, and k = kr + iki is the complex wavenum- 
ber. It is well known that the plane Couette flow is viscously unstable for Reynolds 
number Re > Rcc from the exp eriments but viscously stable from the Orr-Sommerfeld 



equation (jCriminale et aLl 120031 ). If applying the transform in Ea. (l4.2p all the inviscid 
flows must be unstable. Thus the plane Couette flow might be stable only for Re < Rec 
due to the stabilization of the viscosity. 

It is argued that the Taylor- Goldstein e quation r epresents temporal instability, the 
transform represents spatial instability (Hu erre fc Ros si 1998; ICriminale et a/.l2003). I n 



that the pertu rbation is seen along with the flow at the speed of {U — c) in iMilesI (|l96ll ) 



Howard! (|1961[ ). The transform F = (f)/ {U—c)" also turns real wavenumber k into complex 
number, q = implies LOi 7^ 0. The assumption of real k after transform will leads 
to contradictions with the results derived from the Taylor- Goldstein equation. So the 
previous investigators can hardly generalize their results from homogeneous fluids to 
stratified fiuids. 
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5. Conclusion 

In summary, the stably stratification is a destabilization mechanism, and the flow 
instability is due to the competition of the kinetic energy with the potential energy. 
Globally, the flow is always unstable when the total Froudc number Frf ^ 1, where the 
larger potential energy might transfer to the kinetic energy after being disturbed. Lo- 
cally, the flow is unstable as the gradient Richardson number Ri > 1/4. The approach is 
very straightforward and can be used for similar analysis. In the inviscid stratified flow, 
the unstable perturbation must be long-wave scale. This result extends the Rayleigh's, 
Fj0rtoft's, Sun's and Arnol'd's criteria for the inviscid homogenous fluid, but contra- 
dicts the well-known Miles and Howard theorems. It is argued here that the transform 
F = (j)/{U — c)" is not suitable for temporal stability problem, and that it will leads to 
contradictions with the results derived from Taylor-Goldstein equation. 
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Wang W. at OUC and Prof. Huang R-X at WHOI for their encouragements. This work 
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